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Abstract 

The canonical function method (CFM) is a powerful means for solving the Radial Schrdinger 
Equation. The mathematical difficulty of the RSE lies in the fact it is a singular boundary value 
problem. The CFM turns it into a regular initial value problem and allows the full determination 
of the spectrum of the Schrdinger operator without calculating the eigenfunctions. Following the 
parametrisation suggested by Klapisch and Green, Sellin and Zachor we develop a CFM to optimise 
the potential parameters in order to reproduce the experimental Quantum Defect results for various 
Rydberg series of He, Ne and Ar as evaluated from Moore's data. 

PACS numbers: 03.65.-w,31.15.Gy,33.20.Tp 



1 



INTRODUCTION 



A reliable calculation of the electron impact ionisation cross-section of an atom requires 
an accurate determination of the continuum wavefunctions in the incident and in the exit 
channels. At intermediate and low-energy it is necessary to have a fairly good representation 
of the interaction between the incoming electron and the target in the incident channel 
and between the outgoing electron and the residual ion in the exit channel. This work is 
concerned with the latter problem. 

The main difficulty in the calculation of the final state continuum wavefunction arises 
from the non-local character of the exchange potential. Various alternatives have been used 
to overcome this difficulty, one of the most popular being the local exchange approximation 
of Furness and Mc Carthy [jlj. While this approximation seems to work well when the 
continuum energy is high enough it breaks down at low energy. This is precisely the 
case where our model potentials are expected to be the most accurate. 

Our goal is to determine potentials which incorporate the effect of exchange while keeping 
a local character. The potential felt by an electron in the field of an ion is not exactly 
Coulombic at intermediate distances. The difference between the effective and Coulomb 
potentials manifests itself through the Quantum Defects (QD) associated with bound states 
and the phaseshifts associated with continuum states. For a given angular momentum 1, the 
QD oi is defined by: 



where E n i is the binding energy of the excited electron in the nl orbital and Z the 
charge of the ionic core. As n increases the QD is constant for a given Rydberg series. An 
interesting account of the early and recent development of this concept has been given by 



there is in principle a different value of the QD for each coupling scheme between the core 
and the excited electron spin and angular momenta. 

When the electron energy is above the continuum limit, its wavefunction is characterised 



Z 2 



(1) 




Rau and Inokuti [Q. It should be noted that when the ionic core is an open shell system, 



2 



by a phaseshift for each spin and angular momentum. The difference between the actual and 
the Coulombic phaseshifts 5i is the signature of the non-Coulombic part of the electron-ion 
interaction. Near the continuum limit it is related to the QD through: 

Si = 7T<7, (2) 

Thus, our approach is to determine a parametric potential from the information contained 
in the bound state spectrum of the electron-ion system. That potential might be used to 
calculate the continuum wavef unctions. 

An extensive review of the applications of model potentials has been given by Hibbert 
H. Two particularly interesting parametric potentials will be considered here. The first one 
is due to Klapisch ||: 

V(r) = — (2/r)[l + (Z — l)exp{—a,\r) + Cr exp(— a 2 r)] (3) 

where oti, a 2 and C are adjustable parameters. This type of potential has been widely 
used by Aymar et al. (see Aymar || and references therein). The second has been suggested 
by Green, Sellin and Zachor (GSZ) |7| 

V(r) = -(2/r) \{Z - l)O(r) + 1] (4) 

where: 

fi(r) = l/(ei[ezp(-) - 1] + 1) (5) 

where e\ and e 2 are adjustable parameters. Both expressions are Coulombic near the 
origin and at large distances. 

The canonical function method is a powerful means for calculating the QD for a given 
potential since it provides the eigenspectrum of the Schrdinger operator accurately and 
quickly without having to determine the eigenfunctions. 
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THE CANONICAL FUNCTION METHOD 



The canonical function method (CFM) || is a powerful means for solving the Radial 
Schrdinger Equation (RSE). The mathematical difficulty of the RSE lies in the fact it is a 
singular boundary value problem. The CFM turns it into a regular initial value problem 
and allows the full determination of the spectrum of the Schrdinger operator bypassing the 
evaluation of the eigenfunctions. 

The Canonical Function Method (CFM) developed by Kobeissi and his coworkers || [TI| 
to integrate the RSE consists of writing the general solution as a function of the radial 
distance r in terms of two basis functions a(r) and /3(r). Picking an arbitrary point r at 
which a well defined set of initial conditions are chosen ie: a(ro) = 1 with a'(ro) = and 
P( r o) — with /3'( r o) — 1; the RSE is solved by progressing towards the origin and towards 
Infinity During the integration, the ratio of the functions is monitored until saturation 
signaling the stability of the eigenvalue spectrum. 



Several advantages to the CFM are worth mentioning. First, the RSE integration 
belonging to the Singular Boundary Value class is avoided by transforming the problem into 
an initial value one with well-defined starting values for the basis functions a(r) and /?(r). 
Second, the RSE eigenvalue problem is solved and the eigenfunctions are simultaneously 
determined once the roots of the eigenvalue function (EF) are found. The EF is defined as 
the difference of the ratio at saturation for large and small r ||. The EF, depicted in 
Figure 1 below, behaves typically as a tan function when the spectrum has stabilised after 
the saturation of the basis function ratio with r. 



The speed and accuracy of the method have been tested and compared to standard 
integration algorithms in a variety of cases and for a wide of range of potentials. In addition 
the method has been generalised to the phase shift estimation as well to the continuum 
coupled channel problem M [1C 



In potential optimisation problems the method can only be of great value because of 
its speed, accuracy and ease of programming. Accordingly, one has to pick an appropriate 
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optimisation method and use the CFM as a subroutine for a given parametric potential. 



Eigenvalue function versus energy 
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FIG. 1: Behavior of the eigenvalue function with energy. The vertical lines indicate the eigenvalue 
position. 

We choose the well documented Klapisch || and Green et al. @ parametric potentials 
and Helium, Neon and Argon as tesbeds. The experimental values that we start from are 
the levels and quantum defects constituting a stringent test for the validity of the optimised 
potentials. 



RESULTS 

In table 1, we display the Helium (L=0) singlet and triplet results in the Klapisch (entry 
with cti, a?2 and C) and GSZ (entry with e\ and e<i) cases along with the corresponding QD. 
The CFM result is indicated with QD whereas the experimental value is indicated by Exp. 

For Argon states (1=0, 1=1 and 1=2) the optimisation based on the CFM gives with 
Klapisch and GSZ potentials are displayed in table 2. 

For Neon (S=0 with 1=0, 1=1 and 1=2) states using the GSZ potential, we find with the 
CFM the results displayed in table 3. 
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DISCUSSION 



From the above tables a remarkable agreement between the experimental and optimised 
QD is observed in most cases. While a slight agreement might be observed in some instances, 
one should not attribute it to the CFM or to the optimisation procedure but on the large 
uncertainties or limitations observed in the experimental values. The tabulated QD are 
averaged over the full spectrum that might contain 20 or 30 levels obtained from the CFM, 
while the Moore data did not provide more than 10 levels in the same case. Judging from the 
agreement with the level and QD values in some cases, the Klapisch potential description 
might be superior to the GSZ, however it does poorly in the Neon case. We believe the 
GSZ potential yield with our formulation based on the CFM, results that appear more 
adequate to Neon than the Klapisch potential. Nevertheless we cannot rule out entirely the 
Klapisch description simply because its optimisation parameter space is larger and perhaps 
our optimisation procedure did not sample uniformly the entire space. Our work related to 
optimisation of Klapisch potential for Neon is still in progress. 

Acknowledgements: 

This paper is dedicated to the memory of Dr. Hafez Kobeissi (1936-1998). Helpful cor- 
respondance, regarding mathematical aspects of the RSE, with Dr. Jeff Cash are gratefully 
acknowledged. Dr. Stphane Mazevet provided preliminary results for the triple differential 
cross section versus slow electron angle in e-2e experiments etablishing further the validity 
of our optimised potentials. 



[1] Furness J. and Mc Carthy I.: J. Phys. B6 (1973) 2280. 

[2] Rash J.: PhD Thesis, University of Cambridge UK (1996). 

[3] Rau A.R.P. and Inokuti M.: Am. J. Phys. 65 (1997) 221. 

[4] Hibbert A.: Adv. At. Mol. Phys. 18 309 (1982). 

[5] Klapisch M.: CPC 2, (1971) 239. 

[6] Aymar M., Greene C.H. and Luc-Koenig E.: Rev. Mod. Phys. 68 (1996) 1015. 

[7] Green A.E.S., Sellin D.L. and Zachor A.S: Phys. Rev. 184, (1969) 1. 

[8] Kobeissi H.: J. Phys. B15, (1982) 693. 

[9] Kobeissi H. and Fakhreddine K.:, J. Physique II (France), Vol.1, (1991) 38. 

6 



[10] Kobeissi H., Fakhreddine K. and Kobeissi M.: International Journal of Quantum Chemistry, 

Vol XL, (1990) 11. 
[11] Moore C.E.: Atomic Energy Levels, NBS publications (1949). 



7 



TABLES 



s=o 


«i= 6.32313 
a 2 =4.65907 
C= 5.69105 


QD= 0.14 
(Exp)= 0.13 


ei=0.556962 
e 2 = 0.529966 


QD=0.21 
(Exp)= 0.13 


S=l 


ai= 2.00428 
a 2 = 0.998444 
C= 0.654791 


QD= 0.30 
(Exp)= 0.30 


ei= 0.411661 
e 2 =0.688626 


QD=0.29 
(Exp)= 0.30 



TABLE I: Parametric potentials for Helium (Klapisch and GSZ for 1=0, S=0 and S=l) 



1=0 


«i= 9.87780E-02 
a 2 = 1.00056E-01 
C= 7.90032E-02 


QD= 2.51 
(Exp)= 2.16 


£i= 0.115144 
£ 2 = 0.105113 


QD= 2.16 
(Exp)= 2.16 


1=1 


oti= 1.07863 
a 2 = 1.02077 
C= 1.01160 


QD= 1.64 
(Exp)= 1.73 


£i= 0.999785 
£ 2 = 0.981995 


QD=1.64 
(Exp)= 1.74 


1=2 


«i= 0.856959 
a 2 = 1.08286 
C= 10.39858 


QD= 0.30 
(Exp)= 0.36 


£i= 0.117336 
£ 2 = 0.499883 


QD=0.32 
(Exp)= 0.36 



TABLE II: Parametric potentials for Argon (Klapisch and GSZ for 1=0, 1=1 and 1=2) 
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1=0 


ei= 28.2369 
e 2 = 15.7071 


QD=1.32 
(Exp)= 1.32 


1=1 


ei= 5.54229 
e 2 = 1.47477 


QD=0.88 
(Exp)= 0.89 


1=2 


ei= 2.96440 
e 2 = 0.873570 


QD=0.046 
(Exp)= 0.032 



TABLE III: Parametric potentials for Neon (GSZ for 1=0, 1=1 and 1=2) 
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